Semi-central In+In collisions and Brown-Rho scaling 
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In connection with the claim made at the Quark Matter 2005 Conference that the Brown-Rho 
scaling is ruled out by NA60 data we consider dimuon production from semi-central In-|-In collisions 
in a full dynamical model. It is shown that if only a modification of the density-dependent p-mass is 
allowed, the maximum of dimuon invariant mass spectra is only slightly below experimental one. The 
additional inclusion of the temperature-dependent modification factor shifts the spectrum maximum 
toward lower invariant masses making calculation results incompatible with data. A further study 
is needed to disentangle the BR dropping p mass scaling and strong broadening. 



Because of weak interaction, dileptons play an excep- 
tional role among different probes providing information 
on a state of highly compressed and hot nuclear mat- 
ter formed in relativistic heavy ion collisions. Generally, 
a dilepton yield depends on both global properties of 
matter constituents (hadrons and/or quarks, gluons) de- 
fined by the equation of state and individual constituent 
properties related to their in-medium modification. The 
analysis of the e^e~ invariant mass spectra from central 
Pb+Au colhsions at the bombarding energy Eiab—i58A 
GeV, measured by the CERES Collaboration, for cer- 
tain shows an excess radiation in the range of invariant 
dilepton masses 0.2 < M < 0.7 GeV. A possible inter- 
pretation of this excess was given in terms of a strongin- 
medium /O-meson modification (see review articles 0,01 )■ 
Various scenarios of hadron modification were proposed. 
However, low statistics, insufficient mass resolution, and 
large signal/background ratio in the CERES experiments 
did not allow one to discriminate these scenarios, in par- 
ticular those based on the Brown-Rho (BR) scaling hy- 
pothesis j3| assuming a dropping p mass and on a strong 
broadening as found in the many-body approach by Rapp 
and Wambach 0, Q . 

The proposed NA60 experiment with the muon detec- 
tor, zero degree calorimeter and the refined target area 
does not seem to suffer from these CERES shortcomings. 
The first NA60 results for pair production in In+In 

collisions have recently been presented at the Quark Mat- 
ter 2005 Conference Q (see also Fig.l in Ref. 01)). From 
comparison of these results with Rapp's theoretical pre- 
dictions it was declared that the "BR scaling is ruled 
out by NA60" dimuon data whereas the many-body ap- 
proach gives a reasonable agreement. This claim imme- 
diately raised an objection by Brown and Rho. In their 
two letters 0, 0] they noted three points that were ne- 
glected or simplified in the calculation of the dropping 
mass scenario which make it having nothing to do with 
that they believe BR scaling 0,0 to be. 

All three key points stated in fw, ^ may be essential in 
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the description of dilepton production but still have not 
been properly taken into account in the works available 
now in the literature. This also concerns this letter : 
We assume the validity of the vector dominance and do 
not discuss "sobar" excitations. We basically concentrate 
on the influence of the temperature and baryon density 
dependence in the form of scaling. Moreover, nuclear 
interaction dynamics is presented in detail to make clear 
which states mainly contribute to the observed dimuon 
yield as well as the model parameters used. 

The dynamics of heavy ion collisions is treated in terms 
of a hybrid model where the initial interaction stage is 
described hy the transport Quark Gluon String Model 
(QGSM) ^ and the subsequent stage is considered as 
an isoentropic expansion. The latter stage is calculated 
within the relativistic 3D hydrodynamics ^3 allowing 
different equations of state. In our work, the mixed phase 
Equation of State (EoS) is applied. This thermo- 
dynamically consistent EoS uses the modified Zimanyi 
mean-field interaction for hadrons and also includes in- 
teraction between hadron and quark-gluon phases, which 
results in a cross-over deconfinement phase transition. In 
addition to Ref. the hard thermal loop term was 

self-consistently added to the interaction of quarks and 
gluons to get the correct asyniptotics at T >> Tc and 
reasonable agreement of the model results with lattice 
QCD calculations at finite temperature T and chemical 
potential hb 0- 

The ratio between entropy S and total participant 
baryon charge Qb, is shown in Fig. 1 for In+In collisions 
at the impact parameter 6 = 4 fm and bombarding en- 
ergy 158A GeV estimated within QGSM. Being calcu- 
lated on a large 3D grid, this ratio is less sensitive to 
particle fluctuation as compared to entropy itself. Small 
values o{Qb at the very beginning of interaction result in 
large values of the S/Qb ratio. It is clearly seen that for 
tkin > 1-3 fm/c this ratio is practically kept constant and 
this stage may be considered as isoentropic expansion. 

To proceed from kinetics to hydrodynamics, we evalu- 
ate conserving components of the energy-momentum ten- 
sor Too, Toi, Tq2, 2o3 and baryon density (the zero compo- 
nent of the baryon current) within QGSM at the moment 
tkin = 1.3 fm/c in every cell on the 3D grid. This state is 
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FIG. 1; (Color online) Temporal dependence of entropy S 
per baryon charge Qb of participants for semi-central In+In 
collision at Eiab =158^ GeV. 
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FIG. 2: The average energy (solid line) and baryon (dashed) 
densities of an expanding fireball formed in In+In collisions. 
Dotted line shows a contribution of quarks and gluons to the 
energy density. 



treated as an initial state for subsequent hydrodynamic 
evolution of a fireball. The time dependence of average 
thermodynamic quantities is presented in Figs. 2, 3. 

As is seen, the average (over the whole excited sys- 
tem) initial energy density and the compression ratio, 
ns/riQ, are not so high, about 3.2 GeV/fni^ and 1.7, re- 
spectively. These values are by a factor of ~ 2 lower 
than those reached in central Pb-|-Au collisions at the 
same energy in the CERES experiments. Both quanti- 
ties fall down in time and exhibit a weak structure re- 
lated to the softest point of the mixed phase EoS [lH . 
The depicted contribution of quark-qluon degrees of free- 
dom to energy density is not large and the muon creation 



from this phase will be neglected below. The tempera- 
ture evolution (Fig. 3) looks quite similarly. The calcu- 
lated critical temperature at the considered finite /is is 
about 160 MeV, so the evolution starts slightly above 
this value. The end point of the temperature evolution 
curve is about 135 MeV and is governed by freeze-out. In 
our model, the freeze-out occurs locally and continuously 
during the whole evolution. The freeze-out condition is 
that the local energy density (including 6 neighbor cells) 
is below a certain value Sfr which was fixed by describ- 
ing pion multiplicity in central Pb-|-Pb collisions at the 
maximal SPS energy [l^. One should note that the ID 
Lorentz invariant Bjorken hydrodynamics predicts 
an essentially faster fall down and, therefore, a much 
shorter hydrodynamic evolution time. 

To find observable dilepton characteristics, one should 
integrate the emission rate over the whole time-space 
X = (<, x) evolution, add the contribution from the freeze- 
out surface (hadron cocktail), and take into account the 
experimental acceptance. To simplify our task, we con- 
sider only the main channel tttt l^l^ . In this case the 
dilepton emission rate is 

xlm Uern{q,T{x), fibix)), (1) 

where the integration is carried out over the whole grid 
and time from t = till the local freeze-out moment. 
Here = = _ Bose distribution function 

is defined as 

fB{qo,T{x)) = {e'^o/Ti.)_^yl (3) 
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FIG. 3: Evolution of the average temperature. The dotted 
line corresponds to the Bjorken regime with ultrarelativistic 
ideal gas EoS. 
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and lepton kinematic factor is 



CiM) =1 + 2 



1 - 4 — ^- 
M2 



(3) 



with the lepton mass mi. If the pole shift is neglected 
(jyipoie _ mp), the imaginary part of the electromagnetic 
current correlation function, integrated over q, for free 
TTTT annihilation reads 



imn 



g2 (^M2 -rnlY + (Imn) 2 



(4) 



The imaginary part of the p-meson self-energy in the one- 
loop approximation is 



Imn 



(5) 



487r M 

where the parameters are defined in vacuum p^ : 

Qp^^ = 6.05; g = 5.03, nip = 770 MeV. (6) 

To simulate the Brown-Rho scaling, we modify the in- 
medium masses in eq.Q),© estimated according to the 
QCD sum rules by Hatsuda and Lee 0, 0| 



TTLp m*p{x) — mp(l — 0.15 ■ ns(a;)/no) 



(7) 



and simultaneously apply the same modification to the 



1 9* 13, 



as it 



vector dominance coupling m-p/g — > 
was suggested in ,6, 7] . 

Rapp's results cited in ref.0] were obtained using a 
different assumption on the in-medium mass Q 
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The temperature dependence in eq.JHJ is motivated by 
the T-dependence of quark condensate. 

Dimuon invariant mass spectra from In+In collisions 
are presented in Fig. 4 for the tttt channel. Only the 
muon rapidity cut 3 < yiah < 4 is taken into account in 
our calculations. The NA60 experimental points demon- 
strate quite a high resolution in dimuon mass; however, 
their absolute normalization is lost. As is seen, if the 
T-dependent dropping mass © is assumed, our curve 
is strongly shifted toward the low-mass region which is 
in qualitative agreement with the Rapp result 0. In 
the case of the density dependent p mass scaling {Tjl the 
maximum position is only slightly below experimental 
ones and this phenomenological scenario cannot unam- 
biguously be " ruled out" . To compare properly dileptons 
from the pion annihilation with experiment, the contribu- 
tion of free p-meson decay should be added. A number of 
these /9-mesons was estimated as thermodynamical emis- 
sion from frozen-out cells, according to the general pro- 
cedure described in In the last case the agreement 
with experiment will be even better if the muon contri- 
bution from the free p decay is taken into account. Note 
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FIG. 4: (Color online) Invariant mass distribution of dimuons 
from semi-central In-|-In collisions at the beam energy 158 
AGeV. Experimental points are from Solid and dashed 
curves are calculated using the p-mass modification factors 
||7|l and JHJ, respectively. The dash-dotted curve neglects any 
in-medium modification. The dotted line indicates the hydro- 
dynamically calculated p- meson decay at the freeze-out. 



that the calculated solid and dotted lines in Fig. 4 are 
multiplied by the same normalization factor and, there- 
fore, the shown ratio between the signal and free p-meson 
decay is kept. Muon emission from pion annihilation in 
the vacuum case very nearly follows the measured points; 
however, as noted above, the observable quantity is not 
this spectrum but its sum with the hadron decay cocktail 
(dotted line in Fig. 4). It is of interest that making use 
of the not modified m^pj g ratio increases the muon yield 
by a factor of ~ 2, as expected by -S, ,7], but it is not 
visible in normalized data because the spectrum shape is 
weakly changed. 

Finally, the in-medium modification of the p-meson 
mass based on the T-dependent scaling ((SJ is hardly com- 
patible with the NA60 data. As noted in Q , a more 
realistic case corresponds to neglecting any temperature 
dependence. This is indeed so, as follows from our re- 
sults. To disentangle BR dropping p mass scaling Q 
and strong broadening i^j^J, further detailed theoretical 
investigations (e.g. see |l9| ') and a comparison with mea- 
sured absolute muon yield taking properly into account 
the experimental acceptance are needed. 
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